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oT Abstract 



A coarse-grained computational model is used to investigate the effect of a fluid membrane on 
patchy-particle assembly into biologically-relevant structures motivated by viral cores and clathrin. 
For cores, we demonstrate a non-monotonic dependence of the promotion of assembly on membrane 
stiffness. If the membrane is significantly deformable, cores are enveloped in buds, although this 
effect is suppressed for very flexible membranes. In the less deformable regime, we observe no 
marked enhancement for cores, even for strong adhesion to the surface. For clarthrin-like particles, 
we again observe the formation of buds, whose morphology depends on membrane-flexibility. 
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In self-assembly, the interactions between a collection of components guide them to spon- 
taneously form an ordered structure Biological self-organization happens within cells, 
from which all living organisms are composed. Cells are all bounded by a membrane, as are 
many sub-cellular structures. Thus many self-assembly processes are membrane-influenced. 
Membranes themselves are also self-assembled, primarily as a lipid bilayer Q. We focus, 
however, on structures assembled only from proteinaeous sub-units, particularly viruses and 
clathrin. 

The genome of a virus is contained in a core or capsid, a typically mono-disperse shell, 
assembled from individual protein complexes. Often the shells are approximately spherical, 
with many having icosahedral symmetry pf. Viruses are divided into enveloped and non- 
enveloped types, depending on whether the core is surrounded by a membrane. The envelope 
in the former group is acquired through budding 4j. For both enveloped fl-Q] and non- 



enveloped llMl3| viruses there is abundant evidence of membrane influence on core assembly. 



Clathrin, on the other hand, is intrinsically linked to membranes: its main function is the 
formation of coated vesicles for intra-cellular protein transport [l^]. Its three-legged shape 
allows a collection of individual units to form structures that range from extended hexagonal 
sheets to closed cages, which always include 12 pentagonal, in addition to different numbers 
of hexagonal, faces 15]. Assembly is nucleated on cellular membranes by adaptors, protein 
complexes which bind the lattice to the membrane. H exag onal sheets on membranes are 
observed [l6| and coated vesicles form through budding 14 1. 



Experimentally, the reversible disassembly and reassembly of viral capsids in solution 
may be triggered by raising and lowering the pH [17J , allowing in vitro experiments of bulk 
assembly, which is observed, for example, by light scattering 18] or electron microscopy frol ]. 



Similar experiments with clathrin 



20( observed bulk assembly into cage structures, finding 



them to be much more homogeneous when adaptor proteins are present. 

Much theoretical work on biological bulk self-assembly has used patchy-particle models. 
Patchy-particles have discrete, attractive interaction sites on their surface and are very ver- 
satile in terms of the range of structures that may be assembled 21|. The main focus has 



been on the assembly of mono-disperse viral capsids [22j-|26|, with simulations reproducing 
key characteristics such as a lag time, hysteresis and partial capsid formation at high con- 
centrations. Simulations also give more detailed, experimentally inaccessible information 



about assembly dynamics. A similar coarse-grained simulation approach was also applied 



to clathrin assembly 27 



Previous applications of coarse-grained models to the effect of membranes on self-assembly 



are limited, although the effect of rigid templates has been considered 
isotropic spherical particles, on fluid 



29(. Aggregation of 



30|, and polymerized [3l|, membranes was studied. 



More detailed models were also applied to self-assembly within a lipid bilayer [32|, |33 |. 
Although not considering assembly, simulations of particles attracted to a membrane also 
saw budding 34J. An alternative continuum approach found that self-assembly induced 
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budding is controlled by interaction strengths and rigidity 

Given the evidence of the influence of membranes on the self-assembly of biological struc- 
tures, it is important to explore the generic physics that plays a role in such systems. 
Although biological detail is undoubtedly important, we choose rather to investigate coarse- 
grained models that share salient features with viral cores and clathrin. We focus on a 
range of interaction strengths that cover the crossover to assembly in the bulk, as well as 
that from unbound to membrane-bound structures. We correspondingly choose bending 
rigidities within a range whose lowest end gives membranes that are easily deformed, and 
whose highest end gives ones that may not be deformed, by the assembled structures. Here 
we focus on equilibrium, postponing dynamics to a later work. 

Our model comprises Ngu assembling sub-units plus the membrane. The former are mod- 



eled as spherical patchy-particles with a Kern-Frenkel potential 



2l| . similarly to previous 



work 25] but modified suitably such that its first derivative is continuous. The membrane is 



represented using a dynamically triangulated surface model [36(: N mem particles connected 
with 3N mem bonds form a network of N tr i = 2N mem triangles. We sample using Monte Carlo 
(MC) simulations j^l, performed in a periodic rectangular box of sides L XJ L y and L z . The 
membrane's projection completely covers the box in the xy-plane, connecting to itself across 
the boundaries. To apply no external tension 
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we allow L x = L y to vary, whilst also 



adjusting L z to keep to volume, V = L x L y L z , fixed. The MC moves used do not allow the 
membrane topology to change. Sub-units interact with the membrane from both sides but 
are only attracted to one side. Quantities are given in units of the thermal energy, ksT, or 
the typical length of a membrane bond, I. 

Interactions between sub-units (ss) and between sub-units and membrane particles (ms) 
are of a Lennard- Jones type (see appendix). An orientational dependence of the attractive 
part creates patches. The parameter 6q defines the maximum angular deviation of the patch 



position from the particle-to-particle vector before the attractive interaction decreases. It is 
chosen such that, for a given pair of sub-units, only one pair of ss patches can interact at 
once. For ss interactions, twisting of sub-units around an interacting patch is also penalized, 



mimicking the torsional constraints in protein-protein interactions [25]. The ss patches 
are evenly spaced around the particle, with the ms patch lying on the axis of rotational 
symmetry. The sub-units have a size of ~ 2.5 and we choose 8q = tt/4 for the ms interactions. 
This gives a relatively wide ms patch, so that a sub-unit typically interacts with many 
membrane particles so seeing a smooth surface. The minima of the ss and ms interactions 
are — e ss and — e ms . A pair of sub-units are defined to be bonded if their interaction energy 
is < — 0.25e ss . The bending stiffness of the membrane is set by A& ~ k, the bending rigidity 



of the membrane (see 



and appendix) 



We choose two different parameters sets. For our core model, Nsu = 12. Sub-units have 



five ss patches with 6 = 0.2, giving a similar patch width to the optimum in Ref. [25]. 
If the twelve sub-units are placed on the vertices of an icosahedron, they may be aligned 
with the ms patches pointing outwards and every ss patch pointing directly at a patch on 
a neighboring sub-unit, bonding with it. Here, V = 1.07 x 10 4 and N mem = 576. In our 
clathrin-like model Nsu — 36 and sub-units have three ss patches each with 9q = 0.3. The 



patches are wider to allow for a range of curvatures. Following Ref. [28j , the ss patches make 
an angle of (79/180)7r to the ms patch so that, if a closed cage is formed, the ms patches 
point inwards. Here, V = 2.15 x 10 4 and N mem = 1156. Membrane sizes were chosen, using 
preliminary runs, to give plenty of area to cover assembled structures, with the V chosen to 
allow large membrane deformations. Qualitative results were not sensitive to V. 

The main connection to biological systems is that the interactions drive our models to 
form similar structures. The core sub-units resemble intermediate capsomers in the assembly 
of a T = 1 capsid, the smallest virus structure with icosahedral symmetry j^]. In reality, 
enveloped viruses are larger. In our clathrin-like model sub-units are considered equivalent 
to one clathrin, with each patch representing a leg. This is a simplification in that, in 
structure formation, multiple legs of different clathrin lie along each other. We neglect 



adaptor proteins 14|. 



Efficient sampling of our system must overcome a number of issues: free-energy barriers 
between assembled and disassembled states; the importance of collective motion for mem- 
brane relaxation; large times to find target structures. MC allows us to combine different ap- 



proaches that address these problems, specifically Aggregate Volume Bias (AVB) moves 39], 
Hybrid Monte Carlo (HMC) moves 40] and Multicanonical Parallel- Tempering (MPT) 41 ]. 
The AVB moves, which shift sub-units directly between non-bonded and bonded states, as 
well as displacing bonded clusters onto or off the membrane, allow target structures to be 
found very quickly. HMC, which uses molecular dynamics (MD) trajectories to create trial 
states, captures collective motion. Finally, the free-energy barrier problem is ameliorated 
through the use of MPT, involving parallel tempering swaps in two dimensions, e ss and 
e ms . The further addition of a one- dimensional biasing potential, w(E ss /e ss ), constructed 
iteratively during initialization, increases the swap acceptance rate. E ss is the total interac- 
tion energy between all sub-units. We found that the HMC acceptance rate is significantly 
improved by constructing w(E ss /e ss ) as a continuous different iable function and including 



the resulting forces in the MD integration 



42]. 
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FIG. 1. Probability of finding a correctly assembled icosahedral core in a simulation with 12 sub- 
units as a function of e ss and e ms for different (a) v3/2 (b) y/3 (c) 8\/3. Results for 2\/3 and 
are shown in the appendix, (d) As a function of e ss and for e ms = 0.84. 



We define a core to be assembled if all twelve sub-units are in a cluster and each makes 
five bonds. In Fig. [TJwe plot the probability of finding an assembled core, P a , as a function 
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of e ss and e ms for a range of X b between y/3/2 and 8\/3. Our chosen e ss range covers the 
crossover from P a « to P a ~ 1 for a bulk system with the same free assembly volume. For 
all Xb, we observe that, for the lowest e ms , this crossover occurs at about the same e ss as in 
the no-membrane system (see appendix). 

For more deformable membranes, as e ms is increased, assembly occurs at lower e ss . This 
enhancement depends non-monotonically on Xb, see Fig. QJd), occurring over a larger area 
of parameter space for = \/3 than for X b = \/3/2, but then reducing and disappearing 
as Xb is increased further. For lower Xb, and high e ms , the membrane tends to envelop the 
sub-units. In Fig. [21 typical configurations observed for Xb = V3/2 and A& = \/3 with an 
assembled core attached to the membrane are shown. Interestingly, whilst for Xb = VS 
this envelopment is almost complete, forming a bud, for Xb = y/3/2 it is only partial. In 
Figs. [2](c) and (d) we plot the average of the total membrane-sub-unit interaction energy, 
(E ms ), as a function of e ss and e ms for the same Xb, confirming that for Xb = VS/2 the 
membrane envelops the sub-units less: for Xb = v3 the minimum of (E ms ) is ~ —100, whilst 
for Xb = V3/2 it is ~ —70. The lowest (E ms ) are strongly correlated with envelopment in 
buds. 

For X b = 2y/3, some configurations with similar envelopment as for Xb = v3 are seen but, 
for higher Xb only some deformation, not full envelopment, is seen (see appendix). The lack 
of an enhancement of assembly in this regime, despite strong attractions to the membrane, 
is in contrast to the case of extended crystals, where structures grow near attractive walls 



even if the bulk is fluid 



43] 



In the bulk, the probability of assembly is determined by whether the attractions are 
sufficient to overcome the associated entropy loss. The attraction of sub-units to the mem- 
brane confines them, reducing this entropy loss, which may promote assembly. If the free 
energy gain in forming the core is sufficient to overcome the bending energy, as well as the 
entropic cost of binding to the core, the assembled structure may form a bud. Budding is not 
necessary for assembly promotion but the membrane stiffness with the most budding also 
has the most promotion. The non-monotonic rigidity dependence may be due to budding 
suppression: for low stiffness by membrane entropy and for high stiffness by bending energy. 
For our cores, changing the membrane stiffness and attraction changes the probability of 
forming one specific structure. For isotropic particles, in contrast, altering these parameters 



may lead to qualitatively different structures 
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Whilst our range of bending rigidities approximately overlaps with that expected for 
biological membranes (~ 2.5 — 25 [44j), those where we see budding are somewhat on the 
lower side of this range (At, < 2^/3). This discrepancy may well arise from the coarse-grained 
nature of our model and particularly from the relatively small number of sub-units in our 
cores: assuming the free energy gain is proportional to the number of sub-units forming 
them, smaller structures will be less able to deform the membrane into a given shape. 




(a) (b) 




E ms [ k B T l E ms [ k B T l 

(c) (d) 

FIG. 2. Typical configurations for an assembled core strongly attracted to the membrane, taken 
from simulations with e ss = 5.55 and e ms = 1 for \: (a) \^3/2 (b) \/3- Membrane particles are 
shown in green and sub-units in yellow. The positions (not extents) of the ss patches are shown in 
red and the ms patches in blue. The average total membrane-sub-unit interaction energy, {E ms ), 
as a function of e ss and e ms is plotted for the same Aj,: (c) y/3/2 (d) V3. 



For the clathrin-like model, the structures formed are typically poly-disperse, see Fig. [3 
and we use a standard measure of asphericity, A, to investigate their shape. A (see 45] 
and appendix) takes values between and 1, with 1 corresponding to a shape with spherical 
symmetry and corresponding to a non-spherical, oblate or prolate shape. 

We first focus on the average of the number of sub-units in the largest bonded cluster, 
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FIG. 3. Typical configuration from simulations with the clathrin-like model: (a) e ss = 12, without 
membrane, (b) e ss = 12, e ms = 1, A fe = \/3 (c) e ss = 12, e ms = 1, A& = 2\/3. (d) e ss = 12, e ms = 1, 
A;, = 8a/3- Coloring as in Fig. [2j 



(N max ). We consider the same range of e ms as for the core model and choose the e ss range 
so that for the no-membrane system it covers the crossover from small clusters of a few 
sub-units to most of the 36 sub-units being in one cluster (see appendix). For higher e ss , 
without a membrane, the sub-units are observed to form closed cages, see Fig. [3] (a). We 
observe that the sub-unit bonds form 5 and 6 member closed rings on the cage surfaces 
but the shape of the "faces" they enclose deviate significantly from pentagons or hexagons, 
being not generally fiat, and the cages, whil st q ualitatively similar, are not generally of the 



form of the structures observed for clathrin 



151 ] . The key difference may be that when two 



clathrin bond their legs lie along each other, which will result in a greater flexibility to tilt 
the symmetry axes of the two clathrin than to rotate around the symmetry axes. In our 
clathrin-like model, however, a bond has equal flexibility for both such deformations. 

In Fig. H] (a) and (b) we plot (N max ) as a function of e ss and e ms for A& = y/3/2 and 
8\/3. Looking at Fig. Ufa) we see that, for the most flexible membrane, there is a similar 
enhancement of assembly for higher e ms as for cores. However, as shown in Fig. H^b), in 
contrast to cores, the enhancement remains as A& is increased. The results for all intermediate 

8 



35 
30 
25 

A 

20 3 

E 

15z 

V 

10 
5 


1 

0.8 

0.6 A 

< 

0.4 v 
0.2 


FIG. 4. Clathrin-like model: average number of sub-units in the largest cluster, {N max ), as a 
function of e ss and e ms for A;,: (a) v3/2 (b) 8\/3. Average asphericity of the largest cluster, (A), 
as a function of e ss and e ms for A&: (c) 2y3 and (d) 8\/3. 

A(, were very similar (see appendix). This is due to the ability of the clathrin-like sub-units 
to form structures with a range of curvatures: as depicted in Fig. [3(b), at lower A& the sub- 
units form roughly spherical structures which enclose membrane buds. Indeed, at A& = VS/2, 
these are often nearly closed, with the membrane in the bud connected to the rest by a very 
narrow neck (see appendix). At intermediate A&, we observe the formation of more open, 
curved structures, or pits, on the membrane, see Fig. [3(c), whereas at high A& the sub-units 
form extended, approximately flat structures lying on the membrane surface, as shown in 

Fig.EKd). 

In Fig. |4](c) and (d) we plot the average of the asphericity of the largest sub-unit cluster, 
(A), for At = 2y / 3 and 8\/3. For low e ms and higher e ss there is a region for all Xb where 
(A) w 0, corresponding to closed cage structures, not attached to the membrane. For higher 
e ms , non-closed structures bound to the membrane are formed. At lower A& these remain 
somewhat spherical, whilst for A& = 8 a/3 they are a lot less so. 

To summarize, we described a simple coarse-grained model to simulate the effect of a 
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membrane on the assembly of proteinaceous sub-units. We used this model to investigate 
the assembly of structures that share key features with viral cores and clathrin. In both 
cases we found that attraction to the membrane may enhance assembly in regions without 
bulk assembly. For cores, this effect shows an interesting non-monotonic dependence on 
membrane rigidity, being reduced for very deformable membranes and disappearing for the 
stiffest, in contrast to extended crystalline structures with attractive walls. For clathrin-like 
particles, the promotion of assembly persists for less deformable membranes. The difference 
to cores is due to the ability of clathrin-like particles to form structures with different 
curvatures. Furthermore, we observed the formation of biologically relevant buds for both 
cores and clathrin-like particles. In the case of cores, we found that these do not occur if 
the membrane is very flexible, whilst for clathrin-like particles their morphology depends on 
membrane- flexibility. 

The formation of buds on membranes is crucial in various biological processes, for exam- 
ple endocytosis, in which, in some organisms, clathrin plays an important role. Endocytosis 
is a complex process involving the collaborative binding of variety of proteins to the mem- 



brane 



461 ] . The demonstration of bud- formation through assembly in our simulations opens 



the possibility that simple, patchy-particle models could capture basic features of such pro- 
he effects described mi ght also be experimentally observec 



cesses, giving new insight, 
by mixing patchy colloids 47], |48] with giant vesicles 49|, |50|, whose bending rigidity 51] 
lies well within the range considered. More generally, our results clearly demonstrate that 
membranes can have a profound impact on self-assembly and will hopefully stimulate fur- 
ther study in this direction. In future it will be interesting to investigate the dynamics of 
membrane-influenced assembly. 

This work was supported by the Austrian Science Fund (FWF): M1367. Snapshots were 



created using VMD [52J. The computational results presented have been achieved in part 



using the Vienna Scientific Cluster (VSC). 
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APPENDIX 



We give additional details of our simulation model. We also include a number of extra 
plots and snapshots. Parameters and variables are as defined in the main text. 

Model Details 

For the interactions between membrane particles we use smooth potentials that are also 
appropriate for molecular dynamics 53J. Bonded membrane particles interact via 



Ubondij'ij j 







for Tij < 1.15/, 



[80k B T) exp[l/(1.15/ - ry)]/(1.33Z - r { 



(1) 



oo 



for 1.15/ < nj < 1.33/, 



for > 1.33/, 



with r. 



ij 



|r, —rj|, where is position of particle i. All pairs of membrane particles 



experience an excluded volume potential 



UEv( r ij] 



OO 



for < 0.67/, 



[S0k B T) exp[l/(r^ - 0.85Z)]/(r y - 0.67/) 



(2) 



for 0.67/ < Tij < 0.85/, 



for 7^ > 0.85/. 



The minimum distance between any two membrane particles is 0.67/ and the maximum 
bond length is 1.33/. 

A unit normal vector is associated with each membrane triangle. Each bond forms the 
side of two different neighboring triangles. Membrane fluidity is included using MC moves 
that attempt to remove a given bond and create a new one between the two vertices of 
its neighboring triangles that were not connected by the original. During this procedure 
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the direction of the normals is always maintained such that if the membrane were in a flat 
configuration all normals would point in the +z-direction. The bending stiffness of the 
membrane is controlled by including a potential i7& enc z = A&(1 — rij • nj) for each bond, where 
rij and rij are the unit normal vectors of the two triangles neighboring the bond and A& is an 
energy. This form for the bending energy has the weakness that the effective bending rigidity 
is shape dependent. We have also performed simulations using an alternative that does not 



share this deficiency 



36]. We observed qualitatively similar results but found simulation 



times were typically increased by about a factor of two. The total membrane area, A, is 
constrained with a harmonic potential, U area = (/cbT)(A — A ) 2 , where A = (a/3 / 4)Z 2 N tri , 
in the Hamiltonian. 

For the case where particles i and j (i ^ j) are both sub-units (ss) or where one is a 
sub-unit and one a membrane particle (ms) the interaction is of the following form, 

Uij = "/area PwCA{rij) + Isidel orient U a tt{rij)} 



k 



for r < r t , 
for r > r t , 



UwcAir) = < 



U, 



for r < r t , 



k 



for r t < r < r s 



att\ 



a(r 



b{r 



(3) 



for r s < r < r c , 



for r > r r 



'26\l/6 



67, 



24192 e 



and b 



387072 e 



The form of 



where r t = 2 1/6 (7, r s = [fr»<r, r c - w s , „ - - Wi? „ - ^. 
U a tt{f) in the range r s < r < r c is a polynomial interpolation used to avoid a jump in the 
potential or its first derivative at the cut-off 54J. The dimensionless factors '"farea, Iside and 
lorient take different forms for ss and ms interactions. For all simulations reported, we set 
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the length scale, a, for ss interactions to a ss = 2.5/ and similarly set a ms = (l/2)(a ss + I). 
The corresponding energy scales, e, denoted e ss and e ms are varied. 

Since the membrane bonds have a relatively broad, flat minimum the membrane particles 
would tend to be locally compressed when an attractive sub-unit is close. For ms interactions 
the 7area factor is used to counter-act this by making the interaction proportional to the 
area that the membrane particle represents: 7 are a = A nei gh/(N ne ighA t ri), where N neig h is the 
total number of triangles that have the membrane particle as a vertex, A ne i g h is their total 
area and A tr i = Ao/NtH- For ms interactions '-/side is used to make only one side of the 
membrane attract sub-units: it takes a value of 1 if the sub-unit is "above" the membrane 
and if it is "below". A sub-unit is determined to be "above" or "below" by finding the 
closest point on the membrane in the ^-direction. If the normal of the triangle enclosing 
the closest point makes an angle of less than ir/2 with the vector from the closest point 
to sub-unit then the sub-unit is "above" the membrane, otherwise it is "below". For ss 
interactions 7 area = j 8ide = 1. 

The centres of the attractive patches are defined by unit vectors. The width of the patches 
is determined by 7 orient, which is a product of functions of the form [55j: 

1 

for < 

COS 2 [(7T/2)(0-0 O )/0l] 



F(0; 0o,0i 



'0, 



(4) 



for 



< 6 < 



+ 



For ss interactions, 7 



orient 



1 orientiy ij 1 ^ii , 



for > 00 + 01- 

= F(0 <; 6 , e 1 )F{9 j] 6 0l 2# , 26,), 



where f2 describes particle orientation. 6i is the angle between the interacting patch on 
particle i and r^-, whilst 8j is between the patch of particle j and — r^. ipij is the angle 
between the projections of the membrane patches of i and j onto the plane perpendicular to 
iij. The factor F(ipij) penalizes the twisting of interacting sub-units. We follow ref. 25j in 
choosing the range for this factor to be double that for the other ones. For ms interactions, 
lorient = lorient^sm, ^s) = F(Q S ; 9o, #1) , where Y sm is the unit vector from the sub-unit to the 
membrane particle, Q s describes the orientation of the sub-unit and 8 S is the angle between 
the membrane patch and r sm . We choose 9, = 0.2 for all interactions. 
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To investigate the shape of assembled structures we employ the asphericity 451 ] : 

A = ^trQ 2 /(trQ) 2 (0 < A < 1), 

where Q is a 3 x 3 tensor whose elements are given by 



(5) 



= ^ Yl - jp Yl r *.« 



' 3, Pi 



(6) 



where the sums are over all the N particles in the structure considered and r i Q is one of the 
three cartesian components of the position of particle i. Q is a traceless counterpart of Q 
defined as Q = Q — AI, where A is the average of the three eigenvalues of Q. 
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Additional Results 




FIG. 5. Probability of finding a correctly assembled icosahedral core, Passemb-, in a simulation with 
12 sub-units: (a) Without membrane as a function of e ss . With the same free assembly volume as 
for the membrane simulations. Error bars show the standard error, calculated from 4 independent 
repeats of the simulation. Lines joining the data points added as a visual aid. With membrane, as 
a function of e ss and e ms , for (b) 2v / 3^_bT' (c) A^/SksT. 
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(c) (d) 

FIG. 6. Typical configurations for the core model for different (a) Configuration showing 
an assembled, enveloped core for e ss = 5.55/cbT, e ms = ksT, Xf, = (b) Configuration 

from the system with the same parameter values as (a), showing a partially-assembled core that 
is also only partially enveloped, (c) Configuration showing some deformation of the membrane for 
e ss = 6.90/c#T, e ms = ksT, Xf, = A^/SksT. (d) Configuration showing no appreciable deformation 
of the membrane for e ss = 6.90/csT, e ms = ksT, Xb = 8\^3ksT. Coloring as in the main text. 
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FIG. 7. Average number of sub-units in the largest bonded cluster, {N max ), for the clathrin-like 
model, (a) Without membrane, as a function of e ss - With the same free assembly volume as for 
the membrane simulations. Error bars show the standard error, calculated from 4 independent 
repeats of the simulation. Lines joining the data points added as a visual aid. (b) As a function 
of e ss and e ms , for = y^/c^T. (c) As a function of e ss and e ms , for = 2\/3ksT. (d) As a 
function of e ss and e ms , for \t> = 4\/3/c£T. 
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FIG. 8. Typical configuration for the clathrin-like model, for e ss = 12ksT, e ms = ksT, Af, = 
\^3/2ksT, showing an almost-closed cage enclosing a region of membrane connected to the main 
part by a narrow neck (on the left of the image, the protrusion on the right is not directly connected 
to the enclosed region). Coloring as in the main text. 
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